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The crystal growth rates of NH 4 H 2 P0 4 , KH 2 P0 4i (NH^SO,, KA1(S0 4 ) 2 -12H 2 0 > NaCl, and gl ■ 

r i?' KC1 NaBr-2H 2 0, (NH 4 ) 2 C1, and (NH 4 ) 2 S0 4 were express, d 
Qmloroofi ‘ h ® fundamental driving force of crystallization calculated from the activity 
l P ri atUr * ated 8olut,0 . ns - ^ he *V netlc parameters were compared with those from the commo 
used kinetic expression based on the concentration difference. From the viewpoii t 
thermodynamics, rate expressions based on the chemical potential difference provide acci r 
kinetic representation over a broad range of supersaturation. The rates estimated usini . 
expression based on the concentration difference coincide with the true rates of crystallize ti 
only in the concentration range of low supersaturation and deviate from the true kineti. s 
the supersaturation increases. 


Introduction 

Crystallization is an important separation and puri- 
fication process used in many areas such as the chemi- 
cal, pharmaceutical, petrochemical, and electronics 
industries. The kinetics of crystallization is essential 
infonnation required for the design of any ciystallization 
equipment. Numerous experimental works investigat- 
ing the kinetics of crystal growth and nucleation have 
been published in the literature. However, a standard- 
ized approach in correlating experimental results to the 
theoretical or semiempirical models to describe the 
kinetic processes has not been established. Conse- 
quently, estimating rates of crystallization from differ- 
ent expressions often leads to inconsistent results. 

The driving force for crystallization is the degree of 
supersaturation which has been commonly expressed, 
for the sake of convenience, as the difference in concen- 
tration between the supersaturated and saturated solu- 
tions. This practice of expressing the crystallization 
rate as a function of the concentration difference causes 
confusion and inconsistency: even dimensionless su- 
persaturations calculated from different concentration 
units result in different numerical values that are not 
proportional to one another, and thus different sets of 
kinetic parameters could be evaluated from a given set 
of experimental data. 

It is well-known that the fundamental driving force 
of crystallization is the difference between the chemical 
potential of the supersaturated solution and that of the 
solid crystal face, which can be used independently of 
units (Mullin and Sohnel, 1977). Using the concentra- 
tion difference in place of the fundamental driving force 
of crystallization is based on the assumption that the 
solute activity coefficient of a supersaturated solution 
can be closely approximated by that of the saturated 
solution, which may cause serious errors in evaluating 
the true kinetics of crystallization. 

Unfortunately, the kinetic expression using the fun- 
damental driving force seldom has been used in crystal- 
lization practices, because there had been virtually no 
experimental data of activity in supersaturated solu- 
tions. Recently, efforts to study thermodynamic proper- 
ties of supersaturated solutions have been carried out 
by a number of researchers utilizing an electrodynamic 
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microparticle levitator (Cohen et al., 1987; Na, 199. 
et al., 1994, 1995). Water activities of many aqu 
supersaturated solutions have been measured usirp 
single micron-sized solution droplet that is electri a_ 
levitated and continuously weighed as the concentr tit 
is increased by a slow evaporation of water. S. It 
activity of the supersaturated solutions can be comp jL 
from the water activity data using the Gibbs -Du if 
relation; hence, we can establish a direct relation d 
between the chemical potential and the concentr* ti 
of a supersaturated solution. 

From this relationship, many kinetic data rept i 
in the literature can now be expressed in terms o- t 
chemical potential difference providing kinetic ex£n 
sions of more exact and thermodynamically accua 
form. The purpose of this paper is to investigate t 
effects of using the fundamental driving force in U 
expression of crystallization kinetics and to quantif t 
errors associated with the conventional use of U 
concentration-based driving force on the actual rat s 
crystallization. 


Theoretical Section 

Driving Force for Crystallization. The fu t 
mental driving force for crystallization can be exp re s 
in dimensionless form: 


(M -fi*)/RT=\n{a/a*) 

For electrolyte solutions, the driving force is expre t 

in terms of the mean ionic activity of the solute, c ; 


V - J<* 
RT 



where c is the concentration and y± is the correspon i 
mean ionic activity coefficient. The fundamenta. t 
mensionless driving force for an electrolyte solution v 
demonstrated to be the same whether the salt dis o 
ates partially or totally in solution (Sohnel and Mu h 
1978). For solutions of hydrate salts, the term ^ 
refers to the chemical potential difference between Li 
hydrate in supersaturated and saturated soluti m 
respectively, even though the existence of a hydrat ■, 
solution may be only hypothetical (Sohnel and Mu ii 
1978). r 
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The use of the dimensionless concentration difference 

o = (c - c*)/c* (3) 

as the driving force in place of A 41 /RT is justified only 
in the case that meets the following conditions: (i) y ± 
* thus ln(a±/a±*) ^ In (c/c*); (ii) a <*: 1, thus ln(a 
+ 1) a, and (iii) v — 1. 

Condition i will be satisfied only for an ideal solution 
or for a solution whose supersaturation is small enough 
so that its concentration is almost the same as that of 
the saturated solution. There are many circumstances 
in crystallization practices in which a high supersatu- 
ration level ( o > 0.1, for example) is encountered and 
thus condition ii is violated. Examples are primary 
nucleation and precipitation processes where relatively 
insoluble materials are produced as a result of reaction. 
Condition iii will be met for only nondissociating salts 
and hence is inapplicable for electrolytes. Therefore, 
using a in place of Aji/RT is an approximation that is 
inadequate in a number of practical situations. 

Solute Activity* Solute activity is computed from 
experimental water activity data as a function of 
concentration. Through the Gibbs— Duhem equation 
relating the chemical potentials of solvent and solute 
at constant temperature and pressure, the water activ- 
ity can be expressed in terms of solute concentration, 
m in molality, for example (Robinson and Stokes, 1959): 
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kinetic expression (eq 7) based on constant parameters, 
k and n, is applicable only in the restricted range of a 
in which the experimental data were obtained to 
determine the parameters and cannot be used for the 
prediction of crystallization rates over a broader range 
of concentration. 

The power law correlation, eq 8, is used here to 
represent the process occurring at the crystal surface 
only. Therefore, the validity of the equation is limited 
by the assumption that the solute concentration at the 
crystal surface is equal to the bulk concentration 
especially at high supersaturation where the effect of 
the mass transfer diffusion resistance is relatively large. 

b. Boundary Layer Diffusion Model. The overali 
rate of crystal grQwth determined both by the mass- 
transfer rate of solute diffusion and by surface integra- 
tion to the crystal lattice has been described in terms 
of the concentration-based driving forces: 



where C] is the interface concentration and can be 
eliminated to give 


R a 



Rtf 


(10) 


-(1 000/MW w ) d(ln a w ) — vm d(ln (y ± m)) (4) 

upon integration of which we get the relation: 

ln(y±/y±*) = 0 - 4>* + f™. dm (5) 
where <p, the osmotic coefficient, is defined as 

. 1000 ^ a w 

0 MW w vm 

Therefore, 

lnCa^a,*) = ln(m/m*) + <fr - <p* + f m $ — 1 dm (6) 

Jm m 

Crystallization Kinetics, a. Power Law Expres- 
sion. For correlation purposes, a simple semiempirical 
power law equation of the form 


which is often numerically evaluated. Equations 9 and 
10 are valid only in a narrow range of concentration. 
Thermodynamically exact overall growth rate can be 
written as: 


or 


-4 '-iff <n> 


R t = k d 'v ln(a ± /a ±i ) = k x 'v ln(a ± /a i *)' 1 ' (12) 


where m is the chemical potential at the interface 
Analogously to eq 10, the overall growth rate, R g , can 
be evaluated numerically by eliminating the interface 
solute activity, a ±l . This overall growth rate, R g , 
expressed as a function of the chemical potential dif- 
ference, presents the true kinetics of growth over a wide 
range of concentration. 


G = kef (7) 

has been used frequently for the expression of surface 
integration rate of crystal growth and nucleation rate. 
The kinetic expression in terms of the fundamental 
dimensionless driving force, as shown in eq 2, is: 


G 




The kinetic orders n and n' would be the same only 
in the concentration range of very low supersaturations, 
where y±/y±* & 1 and ln(a + 1) can be approximated as 
o. The difference between experimentally obtained n 
and n* wijl be greater for the data of a kinetic experi- 
ment conducted at a higher supersaturation range. 
Thermodynamically, the surface integration rate of 
crystal growth is determined by the constant power of 
the fundamental driving force; therefore, true kinetics 
of crystallization can be described by eq 8 only. The 


Results and Discussion 

Kinetics of Crystal Growth, a. Surface Integra- 
tion Kinetics. Experimental data of the crystal growth 
rates of NH 4 H 2 P0 4 , KH 2 P0 4j (NH^SC^, KA1(S0 4 ) 2 - 
12H 2 0, and NaCl at 25 °C and of glycine at 20 °C are 
reported in the literature (Mullin and Amatavivadhana, 
1967a,b; Mullin et al., 1970; Mullin and Garside, 1967; 
Rumford and Bain, 1960; Li and Rodriguez-Hornedo, 
1992). These authors had conducted the kinetic experi- 
ments using fairly high solution velocity in order to 
eliminate diffusional mass-transfer resistance and thus 
to measure the surface integration rates. First, the 
supersaturation data expressed in various units in 
different papers were converted to be in a consistent 
unit of molality (mol of solute/kg of water or mol of 
hydrate/kg of “free” water). The data of growth rates 
were converted to contain the unit of kg of crystal/(m 2 
of surface area*s). 

The water activity data of NaCl, (NH 4 ) 2 S0 4 , and 
glycine solutions reported by Cohen et al. (1987) and 
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Table 1. Coefficients of a.(m) from the Polynomial a w = o 0 + aim + aim 1 4- ...» at 20-25 °C 


ao 

a 1 


02 

a3 


a 4 06 

06 


m range 
of validit 

1.0726 

-3.081 

X 

10- 2 

4.828 x 

lO" 4 -4.008 X 

io- # 

1.629 x 10' 8 -2.627 x 10 11 

0 


3.4-60 

0.9777 

-1.762 

X 

10~ 2 

-1.387 x 

10- 3 9.390 X 

io-‘ 

-2.128 x lO' 8 1.669 x 10‘ 8 

0 


1.6-20 

1.0400 

-1.702 

X 

10" 1 

2.788 x 

10 3 -2.443 x 

lO' 3 

1.055 x 10- 4 -1.767 x 10-« 

0 


0.2-20 

0.9968 

-2.969 

X 

10~ 2 

1.753 x 

10 s -3.253 x 

lO' 4 

3.671 x 10‘ 8 -9.787 x 10‘ 7 

0 


0.1-18 

1.0084 

-4.939 

X 

10~ 2 

8.888 x 

10" 3 -2.167 x 

10' 3 

1.617 x 10" 4 -1.990 x 10 8 

X 

e* 
»-* 
»— « 

1 

10- 7 

0.5-14 

1.0089 

-3.081 

X 

10* 2 

4.828 x 

10" 4 -4.008 x 

lO' 8 

1.629 x 10' 8 -2.627 x lO’ 11 

0 


3.3-25 

1.0008 

-3.531 

X 

10' 2 

2.490 x 

10" 3 -6.729 x 

10" 4 

5.318 x 10" s -8.040 x 10~ 7 

-2.866 x 

10- # 

0.1-15 

0.9975 

-2.173 

X 

10- 2 

-1.053 x 

10' 2 4.253 x 

lO" 3 

-7.780 x lO' 4 6.203 x lO 5 

-1.764 x 

10- 8 

0.1-13 

0.9968 

-2.611 

X 

10~ 2 

-1.599 X l<r 3 1.365 X lO' 4 

-2.317 x lO’ 8 —1.113 x 10" 8 

0 


0.1-23 

0.9996 

-3.116 

X 

10” 2 

-2.112 x 

10' 3 -9.347 x 

10' e 

2.000 x 10‘ 8 -6.47 x 10' 7 

0 


0.1-2.C 


salt 

NH 4 H 2 PO 4 
KH 2 P0 4 

KAl(S04) 2 -12H 2 O fe 
(NH 4 ) 2 S04 
NaCl 
glycine 
KBr 
KC1 
NH 4 CI 
NaBr-2H 2 0 

0 Selected data from Cohen et al., 1987; Na, 1994; and Chiou, 1994. b The molality m is based on the hydrate mol/kg of free waU 
c The molality m is baaed on the anhydrous mol/kg of total water. 



Na (1993) were used to calculate the solute activity as 
a function of concentration. Additional measurements 
of the water activity for supersaturated NH 4 H 2 P0 4| 
KH 2 P0 4 , and KA1(S0 4 ) 2 *12H 2 0 solutions were carried 
out in our laboratory using the spherical void electro* 
dynamic levitator trap. Details of the experiments can 
be found elsewhere (Na, 1993). The water activity 
expressed as a polynomial function of molality (listed 
in Table 1) was used in the numerical calculation of eq 
5 to acquire the activity coefficient ratio as a function 
of concentration. 

The calculated ratio, y±/y±* is plotted as a function 
of dimensionless supersaturation, o (molality-based), in 
Figure 1. The ratio remains very close to unity up to 
the supersaturation degree of o ~ 0.01; as o increases, 
the activity coefficient deviates either positively (for 
NaCl, KA1(S0 4 ) 2 -12H 2 0, glycine) or negatively (for 
NH 4 H 2 PO 4 , KH 2 P0 4 , (NH 4 ) 2 S0 4 ) from that of the satu- 
rated solution. The biggest deviation of the activity 
coefficient is found for the solutions of NaCl and KA1- 
(S0 4 ) 2 -12H 2 0 and the least deviation for (NH 4 ) 2 S0 4 . 

Figure 2 shows the calculated fundamental driving 
force, t^JRT (=v ln(a±/a±*)), as a function of a. The 
big differences between A ^JRT and o for KA1(S0 4 ) 2 *- 
12H 2 0 and NaCl are due to relatively large deviations 
of y± from y±* and the values of v (v = 4 for KA1- 
(S0 4 ) 2 -12H 2 0, v = 2 for NaCl). For (NH 4 ) 2 S0 4 which 
had a relatively small deviation of y± / y±* in all ranges 
of a, the difference is mainly due to the value of v = 3. 
Although it is not apparent from the plots, all the curves 
have a linear section with slope v in the region of a below 
0.01. The leveling down curvature at high o (above 0.5) 
is due mainly to the fact that ln(o 4- 1) becomes 
relatively smaller than a as a increases. 



Figure 2. Fundamental driving forces calculated as function} 
dimensionless supersaturation (molality-based). 

Parts a— d of Figure 3 show the experimental d 1 ,i 
points of the growth rate versus o and versus A jJR T 
The kinetic parameters for the rate expressions of t u 
forms in eqs 7 and 8 were determined from lint a, 
regressions of the logarithmic plots. The evaluate* 
kinetic parameters are listed in Table 2. It should b* 
noted that the kinetic parameters k T and n are good o 
only the experimental range of o indicated in the tal It 
whereas k / and n' are valid for all ranges of supersa u 
ration. The values of n and n are close for the d L 
measured at the range of very low o (less than 0.01) i 
which the fundamental driving force and concentrat 0 
driving force remain more or less proportional to e, c 
other. For example, the values of n and n' are aim >s 
identical for (NH4) 2 S0 4 (measure at o < 0.007) and N ( 
(measured at o < 0.004). On the other hand, KH 2 1 t 
(measured at 0.1 < o < 0.2) shows much different n < i. 
n\ 

The solid lines (extrapolated from the growth 1 
data with respect to A ji/RT) shown in Figure 3t 
represent the true kinetics of crystal growth over a 
ranges of the driving force. The true rates of gro /i 
determined from the fundamental driving force v j i 
plotted against o (broken solid lines). These are tl 
growth rates that would be measured in kinetic exj ej 
ments conducted over a wide range of o . The p o 
clearly show that different values of the kinetic par u 
eters k r and n would be obtained from experim* 
conducted at different ranges of supersaturation. Tl 
dotted lines, which were extrapolated from the grc ^ 
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Dimensionless driving force (a or Ay/RT) 



Dimensionless driving force (<t or Ajj/RT) 


t>- Q n .. onving lorce faorAuW] 

^tions of dhnens ion less* d rfvin^f^^ ^jy m bol s ) Ajj/RT ( hoi Igiw^ sy^ ^The ^ ht “i Tv *”* (< ° at 20 ° C as 

Rumford and Bam. 1960; U and Rodrfguez-Hornedo, 1992.) atavivadhana, 1967, Mullin et al., 1970; Mullin and Garside, 1967a,b; 


— 2. Kinetic Parameter, of Cryatal Growth Evaluated from Experimental Data 



(NH^SO, 

10011 aingle crysUl face ® rowt h Mullin et al., 1970 

KAKSO i iih n measurements using microscope 

KAUS0 4 )2-12H 2 0 fluidized-bed crysUjlizer Mullin and Garaide 1967a b 

glycine fluidized-bed crystallizer R umford and Ba !nf 1960 ' 

(010) face "‘mesa? 81 * 1 fa t CC gT0Wth • Li and Rodrfguez-Hornedo. 

luxujiace measurements using microscope 1992 

rate data with respect to a, are included only to 
illustrate i the inadequacy of the rate expression based 
on a (eq 7). 

In general, the order n obtained from the experimen- 
tal data of low a (<0.01) can be considered a close 
approximation of the true kinetic order n' We are 
interested in how n (determined from experiment) 


1.6 X 10~ 2 2.18 1.3 x> 10' 


« n2J - ? 007 1 9 * 10 ‘ 2 1 05 6 6 x 10* 3 1 05 

0 002-0.007 6.9 2 .20 7 3 x 10' 1 2.20 

0-01 0.13 6.3 x 10-3 j 50 2.3 x 10" 4 1 63 

0.001—0.004 2.6 x 10 2 2.74 5.7 2 74 


0.01-0.08 

0.01-0.08 


5.4 x 10- 4 1.58 3 6 x 10‘ 4 1 61 
5.3 x 10 s 1.41 3.7 x 10-“ L44 

deviates from n' as a increases. The true growth rate 
u can be expressed as: 


G = k'(A ft/RT) n ' = k(a) a" (a) 


( 13 ) 


where k’ and n' are constant, but k(o) and n(a) y are 
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Figure 4. Fundamental driving force versus dimensionless 
supersaturation based on molality. 



Figure 6. n(o)/n' determined from the tangent lines of the plots 
in Figure 4 versus o. 

changing functions of a. Taking the logarithm of both 
sides produces: 

log k' + n' log(A^/RT) = log k(o) + n(o) log a (14) 

Equation 14 indicates that d(log(Aa/i?D)/d(log o) is 
equal to the ratio of the kinetic orders, n(o)/n\ There- 
fore, in the logarithmic plots of Ap/RT versus o (shown 
in Figure 4), slopes of the tangent lines at varying o 
should be equal to n(o)/n\ It is easily seen from the 
figure that the plots have the tangent lines with slopes 
very close to unity when a < 0.01, and the slopes 
decrease with increasing o. The ratio n(a)In determined 
from the slopes with varying o is plotted in Figure 5. 
In all of the crystallizing systems considered here, n(a) 
is within 5% below n in the range of a below 0.1. As 
indicated in the figure, the range of o at which n(a)/n' 
< 0.9 can be considered as the concentration region 
where n(o) significantly deviates from n\ 

b. Diffusion-Limited Crystallization. The overall 
growth rates of KA1(S0 4 )2*12H20 measured by Mullin 
and Garside (1968) at 32 °C were investigated in terms 
of the fundamental driving force using the boundary 



Dimensionless driving force (<ror A^RT) 

Figure 6. Overall growth rate of KA 1 (S 04 ) 2 * 12 H 2 0 at 32 °C 
(diffusion-limited growth) as functions of dimensionless driving 
force a or A/r/7?T. (mean crystal size = 750 //m; solution velocity = 
0.05m/s). The data points represent the surface integration rates 
at 32 °C with respect to o (filled circles) or (hollow circles). 

(Experimental data adapted from Mullin and Garside, 1968.) 

layer diffusion model described earlier. Separate ex- 
perimental data (Mullin and Garside, 1967a,b) were 
used to evaluate the kinetic parameters for the surface 
integration rates and the mass-transfer coefficient with 
respect to o and A jt/RT. The fundamental driving force 
Afi/RT was calculated from the activity coefficient 
function determined at 25 °C. Using the activity 
coefficient function at 25 °C for the data at 32 °C can 
be justified by the fact that the activity coefficients of 
many binary mixtures have a very small dependence 
on temperature in a moderate temperature range 
(Prausnitz et al., 1986). 

Figure 6 shows numerically calculated overall growth 
rates using eq 12 plotted against Afi/RT (solid line) and 
o (broken solid line). The experimental data of the 
surface integration rates are also shown in the figure. 
The false growth rates determined from eq 10 were also 
plotted as a function of o for comparison. It shows that 
the rates determined from eq 10 would be misleading 
in the range of o above 0.1. 

Kinetics of Nucleation. The primary nucleation 
rates of KC1, KBr, NH 4 C1, (NH 4 ) 2 S0 4 , and NaBr-2H 2 0 
at 30 °C are reported by Nyv\t et al. (1970). The kinetic 
parameters had been determined by the measurements 
of the metastable zone width. The nucleation param- 
eters determined by this method are known to give only 
“apparent” kinetics because growth of critical nuclei to 
visible size is neglected. Nevertheless, the resulting 
parameters provide an effective nucleation rate not 
much different from the reed kinetics of nucleation in 
many cases and close enough for engineering purposes 
(Nyvlt et al., 1985). Furthermore, since our purpose is 
simply to illustrate the effects of using the fundamental 
driving force on nucleation kinetics, we are presenting 
the analysis using these parameters. 

The nucleation rate can be written as: 

y 1 (is) 

which is valid only in the experimental o range, and 

J = k n '(A^/RT) m (16) 

which is the fundamental nucleation kinetics valid over 
a wider range of concentration. 








Table 3. Kinetic Parameters of Nucleation of Various 
Compounds 

range of a G ~ G = MAw/RTT 1 ' 

compound measured order m order m' 



range of o 

G — A n o m 

compound 

measured 

order m 

KBr 

0.01-0.04 

2.49 

KC1 

0.01-0.02 

5.98 

NaBr-2H a O 

0.03-0.07 

2.11 

NH 4 C1 

(NH 4 ) a S0 4 

0.01-0.02 

4.57 

0.01-0.02 

4.65 




/:>■ 

/ 



Driving Force of Nucleation (cr) 

Figure 7. Nucleation rates of KBr (plots 1 ), KC1 (ploU 2), NaBr- 
2HjO (plots 3), NH 4 C1 (plots 4 ), and (NH 4 ) 2 S0 4 (plots 5) at 30 °C 
as functions of a. (Experimental data adapted from Nyvlt et al 
1970.) 

After converting the supersaturation data to be in the 
unit of molality, the nucleation parameters were re- 
evaluated in terms of o and AjdRT. The nucleation 
rates were expressed in the unit of kg of crystal/(m 3 of 
aolution-s). The comparisons of the orders m and m' 
are presented in Table 3. 

Figure 7 shows the nucleation rates versus the 
dimensionless driving force a. The extrapolated dotted 
lines representing eq 15 seem to coincide with the true 
nucleation rates depicted by eq 16 (broken solid lines) 
in the range of a below 0.1, but the deviation increases 
with increasing a . 

Conclusions 

Kinetics of crystal growth and nucleation of various 
salts were investigated in terms of the fundamental 
driving force of crystallization based on the chemical 
potential difference. The chemical potentials of super- 
saturated solutions were calculated as a function of 
concentration from the activity of supersaturated solu- 
tions. The kinetic expressions baaed on the chemical 
potential difference are thermodynamically exact rep- 
resentations of the true kinetics of crystallization. The 
commonly used kinetic expressions based on the con- 
centration difference are in significant error at the high 
supersaturation range, in which many crystallization 
processes are commonly operated. 
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Nomenclature 

a = activity of solute 

fit — mean ionic activity of solute 


Ind. Eng. Chem. Key., Vol. 35, No. 4, iyy*> 108c 

a w = activity of water in aqueous solution 
c concentration of solute in any unit 
Cj — concentration of solute at the interface between the 
solution and the boundary layer 
G = crystal growth rate (surface integration), kg/m 2 -s 
J - crystal nucleation rate, kg/m 3 *s 
k = crystallization rate constant, kg/m 2 -s 
k' — crystallization rate constant with respect to A ufRT 
kg/m 2 *s 

= mass-transfer diffusion coefficient, kg/m 2 *s 
k n = nucleation rate constant, kg/m 3 *s 
k n ' = nucleation rate constant with respect to A iJRT k e/ 
m 3 -s * w 

k r - surface integration rate constant, kg/m 2 -a 
k t = surface integration rate constant with respect to A ul 
RT , kg/m 2 *s • K ^ 

m = molal concentration of solute, mol/kg of water 
m — order of crystal nucleation rate expression based on a 
m' = order of crystal nucleation rate expression based on 
Afi/RT 

MW w = molecular weight of water, g/mol 
n = order of crystal growth rate expression based on a 
n - order of crystal growth rate expression based on A/V 
RT 

R = gas constant, J/mol*K 

R t = overall rate of diffusion-limited crystal growth kg/ 
m 2 *s * 

T = temperature, K 
4> = osmotic coefficient 
V± = mean ionic activity coefficient 
H = chemical potential, J 

Hi — chemical potential at the interface between the 
solution and the boundary layer, J 
A/* = chemical potential difference between supersaturated 
and saturated solution, p - p* t J 
v = number of moles of ions formed from 1 mol of electrolyte 
a = relative supersaturation, (c - c*Vc* or (m - 
* = of saturated solution 
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